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Within the universal zero-range theory, we compute the three-body recombination rate to deep 
molecular states for two identical bosons resonantly interacting with each other and with a third 
atom of another species, in the absence of weakly bound dimers. The results allow for a quantitative 
understanding of loss resonances at finite temperature and, combined with experimental data, can 
be used for testing the Efimov universality and extracting the corresponding three-body parameters 
in a given system. Curiously, we find that the loss rate can be dramatically enhanced by the resonant 
heavy-heavy interaction, even for large mass ratios where this interaction is practically irrelevant 
for the Efimov scaling factor. This effect is important for analysing the recent loss measurements 
in the Cs-Li mixture. 
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I. INTRODUCTION 

Measuring atomic three-body losses near Feshbach res¬ 
onances for large and negative scattering lengths has be¬ 
come a major tool for characterizing the Efimov physics 
in a large variety of ultracold gases, not only for identi¬ 
cal bosonic atoms [IHIl] but also for homonuclear [la¬ 
in] and heteronuclear mixtures [IHHII]- The peaks in 
the three-body loss rate as a function of magnetic field 
mark passages of Efimov trimers through the free-atom 
scattering threshold. If the size of such a trimer is much 
smaller than the typical de Broglie wavelengths in the 
gas, the corresponding peak is most visible and is well 
described by zero-temperature theory Then, accord¬ 
ing to the Efimov discrete scaling invariance, the next 
peak is expected to occur when the two-body scattering 
length is multiplied by the so-called Efimov period, the 
trimer being proportionally larger. The Efimov period 
is numerically quite big, so that when trying to observe 
multiple successive peaks, i.e., trying to test the Efimov 
scenario, one very soon faces the problem that Efimov 
states become too large and the loss rate saturates to a 
constant value unmiiiiiiiii. 

For identical bosons the Efimov period equals 22.7 [23] 
and two loss peaks separated by approximately this fac¬ 
tor have recently been observed in Cs [HIS]. However, 
the second (excited-state) peak is already close to the 
saturation regime and its quantitative characterization 
has been done relying on the finite-temperature theory 
developed in Ref. m for three identical bosons. Simi¬ 
larly, an Efimov loss feature for large negative scattering 
lengths in the mixture of three hyperflne states of ®Li 
has been recently reanalyzed by using a suitable mod¬ 
ification of this finite-temperature theory, leading to a 
better determination of the three-body parameter in this 
system m- 

The Efimov period is significantly smaller for strongly 
mass-imbalanced heteronuclear systems [24] . making 
them ideal candidates for testing the Efimov discrete scal¬ 


ing invariance [25] . Very recently, experiments 
have observed more than one period of the Efimov scal¬ 
ing dependence in the system of one Li and two Cs atoms, 
where the Efimov period is ~ 5. For ®Li-®^Rb-®^Rb this 
quantity is ~ 7 and the Rb-Li mixture is thus also poten¬ 
tially interesting from the viewpoint of testing the Efimov 
scenario. 

In this article we generalize the ^-matrix formalism 
developed for three bosons (case BBB) in Refs. [TOl HBj- 
[25] to the case of an atom interacting with two identical 
bosons (case ABB). We assume that the AB interaction 
is tuned to the negative side of a Feshbach resonance and 
consider two cases for the BB interaction: i the BB in¬ 
teraction is neglected, and ii the BB scattering length is 
large and negative. Case i is suitable for the Li-Rb-Rb 
system, while case ii is relevant for the Li-Cs-Cs one in 
the magnetic held region studied in Refs. [2ni[2T]. In none 
of these two cases the system supports weakly bound 
dimers and, therefore, the loss is entirely determined by 
the recombination to deeply bound states. Under these 
conditions we express the loss rate constant in terms of 
the temperature, three-body parameter, inelasticity pa¬ 
rameter, and a quantity sn which is a function of a single 
variable kuAB in case i and of two variables, kuAB and 
kuBB, in case ii, where k is the three-body collision mo¬ 
mentum. Once this universal function is calculated (and 
tabulated) for a given AB mass ratio, one can then eas¬ 
ily generate loss curves for any given temperature, three- 
body parameter, and inelasticity parameter. 

The article is organized as follows. In Sec. [n] we intro¬ 
duce basic notations, present a formula for the loss rate 
constant, and apply it to ®Li-®^Rb-®^Rb and ®Li-^^^Cs- 
^^^Cs commenting on the role of the BB interaction. In 
Sec. [In] we give a detailed derivation of this loss rate for¬ 
mula. The calculation of the universal function sn is 
the subject of Sec. [IV] Appendices include the normal¬ 
ization constant and contact parameters for Efimov ABB 
trimers, the analytic expansion of sn near unitarity, and 
in App. 1 we show how one can recover the case BBB 
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from the present article. 


II. MAIN RESULTS 

We start with some reminders and notations about the 
three-body problem of two bosons of mass tob and a third 
particle of mass uia- We write m for twice the reduced 
mass, 


m = 


mAfriB 
niA + rriB 


( 1 ) 


The mass ratio is conveniently parameterized by the an¬ 
gle 


sin(}i = 


wb 

rriA + wb 


For further convenience we also define 


9 = 


TT (j) 
4 ~ 2’ 


( 2 ) 

( 3 ) 


We denote by oab and obb the AB and BB scattering 
lengths. 

Universal properties, such as the Efimov discrete scal¬ 
ing invariance, appear in the zero-range limit, that can be 
described by a universal theory mm- Accordingly, we 
consider zero-range interactions between A and B parti¬ 
cles, while between B particles we assume no interactions 
in case i and zero-range interactions in case ii. The valid¬ 
ity conditions for this universal zero-range theory are the 
following. In case i, the characteristic interaction ranges 
and IobbI should be much smaller than |oab| and the 
typical de Broglie wavelengths. In case ii, the charac¬ 
teristic interaction ranges should be much smaller than 
IoabIj IobbIj and the typical de Broglie wavelengths. 

The unitary limit is defined by 


Note that we do not consider here the case of a very 
large mass ratio, where additional Efimovian sectors 
would appear in non-zero total angular momentum sub¬ 
spaces. We thus assume that wb/toa is smaller than 
the critical value ~ 40 where an additional Efimov effect 
appears in the angular momentum L = 2 subspace |24j . 

The Ehmov effect also implies that the interactions 
within the universal zero-range theory are parameter¬ 
ized not only by the scattering lengths, but also by a 
three-body parameter. For this purpose we use a length 
parameter i?o defined modulo multiplication by 
This parameter fixes a hyperradial node of the three- 
body wavefunction at small interparticle distances and 
thus hxes all other three-body observables [see Eq. (C 6 ) 
for its relation to the spectrum of Efimov trimers]. To 
account for recombination losses to deeply bound states 
in the zero-range theory one allows the three-body pa¬ 
rameter to be complex, i?o —>■ i?o where rj > 0 

is the inelasticity parameter [1]. A precise definition of 
Rn and n is given by the three-body contact condition 
Eq. @ below. 

Summarizing, the external parameters of the zero- 
range theory are the scattering length oabi three-body 
parameter Rq, inelasticity parameter rj, and also the scat¬ 
tering length obb in case ii. Experimentally, all these 
parameters depend on the magnetic field B. However, 
essential for universal Efimov physics is only the reso¬ 
nant enhancement of oab near a Feshbach resonance. In 
its sufficiently narrow vicinity the other parameters can 
be assumed constant since their H-dependence is smooth. 
It is important to note that the assumption of constant 
Rq and rj is not directly related to the applicability condi¬ 
tions of the zero-range approximation mentioned above. 
In particular, all our derivations and formulas remain 
valid for H-dependent i?o and 77 . However, these param¬ 
eters are assumed constant in Figs. 1 and 2. 


OAB = 00 in case 7 , 

OAB = obb = 00 in case if. 


( 4 ) 


As usual, scattering length values -l-oo and —00 are equiv¬ 
alent. At the unitary limit, the scattering length(s) do 
not introduce any lengthscale into the problem. Due to 
the Efimov effect, there is no continuous scale invariance, 
but a discrete one. In particular, the energies of Efimov 
trimers follow a geometric series, 

A„+i/A„ = (5) 

Here Sq is the positive real solution of [Mj 


A. Loss rate constant 

Our system of interest is a mixture of A and B atoms, 
in the gaseous non-degenerate regime where the thermal 
wavelength is small compared to the interparticle dis¬ 
tances. We do not consider hnite positive values of the 
scattering lengths, where weakly bound dimers would ex¬ 
ist; we also assume that oab is non-zero, in order to have 
a non-trivial three-body problem. 

The rate of ABB recombination events per unit vol¬ 
ume is K ua ^, where ua and riB are the atom number 
densities, and K is the event rate constant, for which we 
obtain 


X{so,4)) = 1 in case i, 
1 — A(so, 4’) ~ 2 A^(so, 0 ) = 0 in case ii, 


where 


A(so,<5) = 


2 sinh(^so) 


So cosh(7rso/2) sin(2(5)' 


( 6 ) 

( 7 ) 


( 8 ) 


K = 647r^ cos^cf) 






(1 - e-4’’) 

(9) 


Jq \1 + {kRo)-'^^^oe-^v Sii\^ 

The function Sn depends only on kaAB, kaBB, and the 
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mass ratio. Therefore, Eq. 0 gives an explicit depen¬ 
dence of the loss rate on the three-body and inelasticity 
parameters. This is a manifestation of a general property 
known as Efimov’s radial law mui]. 

For high temperatures or large scattering lengths such 
that in case i we typically have — fcoAB 3> 1 and in case ii 
—fcoAB S> 1 and — fcoBB ^ 1, the function sn in Eq. ([^ 
can be approximated by its asymptotic unitary value 

gTC _ _g-7rso-|-i2[so ln2-|-argr(l+iso)] 


with So given by Eqs. 

In Section |IV| we show how to compute sn at finite 
/coab and fcoBB for any given mixture. Once this task is 
accomplished, Eq. 0 offers a very fast way of calculating 
the loss rate for any T, iioj and 77 , and can thus be used 
for extracting these parameters from experimental data 
in the universal limit. Obviously, in this manner not 
only experimental but also theoretical results obtained 
for finite-range potentials can be compared to the zero- 
range theory. 

In the next two subsections we employ Eq. 0 for cal¬ 
culating K in two experimentally relevant cases charac¬ 
terized by large mass ratios: A = ®Li, B = ®’^Rb and 
A = f^Li, B = i33cs.[^ 


B. Lithium-Rubidium 

We first consider the case A = ®Li, B = ®^Rb and ne¬ 
glect the ®^Rb-®^Rb interaction since it is non-resonant 
apart from narrow magnetic field intervals [Miisn]. Broad 
interspecies Feshbach resonances in the hyperflne ground 
states are available for both ®Li-®^Rb and ^Li-®^Rb com¬ 
binations [5DII32) . However, for the purpose of test¬ 
ing the Efimov scenario in the ABB system the choice 
of fermionic A allows one to neglect concomitant AAB 
losses suppressed due to the fermionic statistics [55] . 

The ®Li-®^Rb-®^Rb system is characterized by sq = 
1.63188 and a scaling factor = 6.85610. In Fig. 
we present the event rate constant K in units of ha'^/m 
as a function of a-ja for different temperatures. Here a 
is the Rb-Li scattering length and a_ < 0 is the position 
of the three-body loss peak at zero temperature (at a = 
a- an Efimov trimer crosses the three-atom threshold). 
The (three-body) parameter a_ is convenient when there 
is only one relevant scattering length (for instance, for 
three identical bosons or for our case i). We find that 
for the ®Li-®^Rb mass ratio a_ is related to Rq by o_ « 

-5.233Ro- 

The shapes of the curves are controlled by rj and 
T. Unaware of experimental results, we set rj = 0.2. 
The lowest curve in Fig. [l] corresponds to T = Tq = 


^ Tabulated values of sn for these systems can be obtained by 
contacting the authors. 



FIG. 1. The event rate constant K in units of ftal/m versus 
a-/a for the ®Li-®^Rb-®^Rb system for various temperatures 
with the inelasticity parameter set to rj — 0 . 2 . 


li^/kBiTia^ and the other ones are obtained by decreas¬ 
ing InT by 27r/3so, be., a third of the Efimov energetic 
period. Thus, the solid curve for T = Tq and the dotted 
one for T = are self-similar, the latter can be 

obtained by shifting the former horizontally by 
and vertically by The zero-temperature result 

(thin solid) is known analytically (55]. The shift of the 
resonance peak towards smaller |a| with increasing T is 
likely to be related to the fact that for |a| < |a_| the 
trimer becomes a finite-energy resonance in the three- 
atom continuum. 


C. Lithium-Cesium 

Because of its larger mass ratio, the combination 
A=®Li, B=^^^Cs is considered as an even better can¬ 
didate for observing several periods of the Efimov dis¬ 
crete scaling. Two groups have recently reported loss 
measurements for this mixture close to a wide CsLi Fes¬ 
hbach resonance at 843G [201 ED- The first two Efimov 
resonances and signatures of the third one have been de¬ 
tected, although the latter is rather strongly thermally 
saturated. 

In treating the ®Li-^^®Cs-^^^Cs system for magnetic 
fields studied in [551ET] one has to have in mind that the 
CsCs scattering length is rather large (about -1550aBohr 
at the CsLi resonance [5S]) and this may increase the loss 
rate constant because of the enhanced probability to find 
two Cs atoms close to each other. On the other hand, 
to describe the Efimov effect for such a large mass ratio 
one is tempted to use the Born-Oppenheimer approxi¬ 
mation [36| . within which the heavy-heavy interaction is 
irrelevant for the discrete scaling. Indeed, the ®Li-^^^Cs- 
^^^Cs system in case i is characterized by a scaling factor 
g 7 r/so _ 4,87661 [Eq. 0 gives sq = 1.98277], whereas for 
aCsCs = 00 this quantity becomes 4.79887 [sq = 2.00308 
as follows from Eq. 0 ] (cf. [37]). 

In Fig. we plot K vs — I/acsLi for T = 400nK (solid 
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FIG. 2. K vs. —1/acsLi at T = 0 and T = 400nK includ¬ 
ing and excluding the CsCs interaction. The three-body and 
inelasticity parameters are set to Ro ~ ISOoBohr and rj = 0.6, 
respectively. 


and dashed) and T = 0 (dash-dotted and dotted) with 
i?o = ISOoBohr and rj = 0.6. Including the CsCs in¬ 
teraction we obtain the solid and dash-dotted curves 
whereas neglecting it gives the dashed and dotted lines. 
In the CsCs-interacting case we use the magnetic field de¬ 
pendence acsCsiB) provided in Suplemental Material of 
Ref. [35] and eliminate the magnetic field by using the for¬ 
mula acsLi = -28.5aBohr[61.4G/(R-842.9G)-hl] HOj- A 
different acsLi(.S) was used in Ref. [HI and Ref. [35] gives 
yet another more recent characterization of acsuiB). For 
the results plotted in Fig. [^ these variations are not very 
important because acsCs{B) is smooth and, therefore, 
the dependence acsCs(acsLi) is practically unchanged. 

Figure [^ clearly shows that the effect of finite T is to 
saturate K at large acsLi whereas the inclusion of the 
CsCs interaction leads to a strong enhancement of the 
loss rate for small acsLi- For the experimentally studied 
region of B either of these effects can lead to corrections 
of two orders of magnitude. The overall behavior of K 
is thus much less steep than the scaling K oc (com¬ 
pare solid and dotted lines in Fig. [^, consistent with the 
experimental findings [ini[ii]. The peak positions move 
when we “switch on” the CsCs interaction, but this is an 
artefact of choosing the same Rq in these two cases. In 
fact, even in the limit acsLi oo the loss features are not 
expected to exactly match. This is because the Efimovian 
type-i* wave function has to propagate through distances 
of order acsCs before it can be matched with the type-i 
one. In general, when both cases i and ii are applicable, 
i.e., when |obb| is much larger than the characteristic 
interaction ranges but smaller than IoabI and 1/fc, the 
corresponding three-body parameters are related to each 
other by the wave-function matching condition, but are 
not necessarily equal. 

The parameters i?o = ISOoBohr and rj = 0.6 used for 
plotting Fig. [^ have been chosen as a result of a very 
approximate fitting of the data of Refs. A sig¬ 

nificantly more serious account of experimental uncer¬ 


tainties and cross-correlations is needed to give a more 
definite answer for Rq and rj. We find that the fitting 
procedure is very sensitive to the exact position of the 
Feshbach resonance. Surely, one would benefit from more 
experimental data at lower temperatures. Setting T — 0, 
r] ^ 1, and Rq = ISOoBohr we obtain sharp peaks of K 
at positions acsLi = « -330aBohr, acsLi = « 

-1.7 X lO^OBohr, acsLi = « - 8.8 x lO^UBohr, which 

compares very well with the refined analysis of the ex¬ 
perimental data performed in Ref. |38j . 

We note that the experiments [20] [H] are well within 
the zero-range limit, at least, for the second and higher 
Efimov resonances. Indeed, out of the van der Waals 
ranges [Ml 133, the Feshbach resonance parameters i?* 
m 14(11 l4T] , and the length R^,\abg\ that enters when 
the background scattering length Ubg is negative [32] , the 
largest is the Cs-Cs van der Waals range rvdw.CsCs = 
lOlttBohr- In turn, this quantity is much smaller than 
the thermal wavelengths (>3400aBohr at T < 800 nK) 
and the Cs-Cs scattering length, which, for experi¬ 
mentally relevant magnetic fields, varies in the interval 
— ISSOttBohr < ocsCs < —OGOttBohr [35] • As far as the Cs- 
Li scattering length is concerned, around the first Efimov 
resonance it approximately equals — 300aBohr, which is 
comparable to Tvdw.CsCs- However, for the second and 
higher resonances the inequality |acsLi| ^ ?'vdw,CsCs is 
well satisfied. 

Wang et al. [43] have studied the so-called van der 
Waals universality of the three-body parameter in het- 
eronuclear systems close to a wide interspecies resonance 
by assuming the Lennard-Jones interatomic potentials. 
According to their analytical estimates in the case of a 
large mass imbalance the three-body parameter equals 
the heavy-heavy van der Waals range rvdw.BB times a 
dimensionless function of the ratio r-vdw.BB/oBB- For 
the Cs-Cs-Li case in the experimentally relevant region 
this function varies very little (< 2 %) since the ra¬ 
tio rvdw.CsCs/acsCs stays small. We mention for ref¬ 
erence that for acsCs = 2000aBohr Hef. [43] predicts 
a_ = —1.4 X lO^OBohr- 

A finite-temperature theoretical analysis of losses in 
this system has been performed by Y. Wang and reported 
in the Supplemental Material of Ref. m- Wang’s re¬ 
sults indicate that with increasing temperature the reso¬ 
nance features become weaker and shift towards smaller 
oCsLi, which is what we also observe. In fact, we can 
rather well fit Wang’s curves for all three considered tem¬ 
peratures (T = lOOnK, 250nK, and IfiK) by choosing 
Rq = lIOoBohr and rj = 0.4. 


III. DERIVATION OF THE LOSS RATE 
CONSTANT 


In this section, we derive the expression ([^ for the loss 
rate constant in terms of the quantity sn, as well as the 
analytical result (101 for sn valid at the unitary limit. 


Let us first introduce the three-body scattering state 
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ij}. Denoting the positions of the two identical bosons by 
ri and ra, and the position of the third particle by r 2 , 
the wavefunction ^(ri,r 2 ,r 3 ) is symmetric with respect 
to exchanging ri and ra. 

Let us introduce the Jacobi coordinates 


• the three-body contact condition: 3B such that 


^(R) 


Ft —^0 




B{n) 

R2 ■ 


(18) 


cos (/) X = ra — 


mAr2+mBri 

rriA+mB ’ 


y = r 2 - ri. 


( 11 ) 


All information about the relative positions of the three 
particles can then be collected in the 6-dimensional vec¬ 
tor 


The loss rate can then be obtained from the scattering 
state thanks to the following exact relation, that can be 
justified by heuristic arguments. Introducing the proba¬ 
bility current 


J = - Im(^* VrV>) 
m 


(19) 


R = (x,y). 


(12) and the lost flux 


Its norm R = ^is the so-called hyperradius, 
while its direction 


^\o 


= -<t 3 -d^S 


( 20 ) 


n = R/R 

can be parameterized by five hyperangles. One can note 
that if all particle coordinates are multiplied by a fac¬ 
tor A, then the hyperradius is multiplied by A while the 
hyperangles are unchanged. 

At large distances, the scattering state asymptotes to 
an incoming plane wave plus a scattered wave. More 
precisely, in the center-of-mass reference frame, 

^(R) cs V'^°^(R)+V'sc(R) (13) 

i?—>oo 

where ipsciR) is a purely outgoing scattered wave, and 


where 5 is a hypersurface enclosing the origin (e.g. a 
hypersphere) and the surface-element vector d®S points 
away from the origin, the loss rate constant is given by 
the thermal average 


/ (fk 

j e-k^/’^ksT ^ 6 ^ 


of the energy-resolved event rate constant 


K{k) 



(‘/^loss) 


( 21 ) 


( 22 ) 


^/>(°)(R) 


l + P 
72 


g*k.R 


(14) 


is a symmetrized plane wave, normalized in a unit vol¬ 
ume. The operator P in Eq. (14) exchanges particles 1 
and 3. The relation between k and the collision energy is 



(15) 


Here and in what follows, we set h = 1. 

Furthermore, if) is an eigenstate of the zero-range 
model, i.e., it satisfies 


• the Schrodinger equation 

- = Ei,{R) 


(16) 


when none of the three particle positions coincide; 


• the two-body contact conditions: for each interact¬ 
ing pair of particles i,j, 3Aij such that 







Aij -I- o(l) 


(17) 


where Aij depends on the relative position of the 
third particle with respect to the center-of-mass of 
particles i and j ; 


where {■)t denotes the average over the direction of k. 
In Eq. (22), the factor 1/2 originates from the indistin- 
guishability of the two B particles, and the factor cos^cj) 
from the Jacobian of the change of variables from Carte¬ 
sian to Jacobi coordinates, 


d(ri,r2,r3) 

9(C,x,y) 


cos^cj) 


(23) 


where C is the center-of-mass of the three particles. 


A. Unitary limit 


Let us first consider the unitary limit Eq. ([^, where 
the situation is particularly clear, because the problem 
can be solved in a fully analytical way. The key ingre¬ 
dient is that there is a separability between the hyper¬ 
radius R and the hyperangles fJ, because the two-body 
contact conditions © do not introduce any lengthscale 
and hence act only on the hyperangles [24] . The solutions 
of the hyperangular part of the three-body problem are 
the functions (/)s(J^) that satisfy the two-body contact 
condition and which are eigenfunctions of the Laplacian 
operator on the hypersphere: 

To 4>s{^) = -s^(/s(r2). 


(24) 
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The operator Tfj is defined as the hyperangular part of 
the total Laplacian 

We can then expand the scattering state as 

= (26) 

S 

Indeed, the functions (j)s, normalized to unity, form an 
orthogonal basis for the hyperangular scalar product 

{f\g) = I finrg{n)dn (27) 


has the meaning of a reflection amplitude from the point 
R = 0 (where all particle positions coincide). While the 
phase of this reflection amplitude is determined by the 
three-body parameter Rq, its modulus is determined by 
the inelasticity parameter g, the reflection probability be¬ 
ing Accordingly, 

^loss = (1 - e-^'') |ArT (33) 


[as obtained by taking a vanishingly small hypersphere 
for S in Eq. @]. 

In what follows we denote by the long¬ 
distance amplitudes [see Eq. (p^j. The 

out-amplitudes can be expressed in terms of the in¬ 
amplitudes through a linear relation, 


(where dfl stands for the differential solid angle in six- 
dimenional space, i.e. d^R = dfl R^dR). This follows 
from the self-adjointness of the zero-range model. 

In the present case, where the B particles are bosonic 
and the mass ratio is not very large, the set {s} contains 
a single imaginary value s = isq where sq solves Eq. 
and an infinite countable set of real values. The s — 
isQ sector is called Efimovian, since it causes the Efimov 
effect. 

The hyperradial Schrodinger equation reads 



f Sll Sl3\ f \ 

VS31 S33J \A* 3 "y ■ 


Combining this with Eq. (311 yields 


A 


out 

1 


Sl3 

I — S 11 .A 


A 


in 

3 • 


(34) 


(35) 


The 5'-matrix is easily computed by using the fact 
that Fisg{R) is a linear combination of the Bessel func¬ 
tions J ±is„ (kR). This yields the expression of sn given 
in Eq. (10). Furthermore the obtained matrix is unitary, 
so that 7 si 3 p = 1 — |siiP and thus 


i.e., the unitary three-body problem reduces to a set of 
independent one-body problems in effective po¬ 

tentials. The hyperradial wavefunctions Fs{R) have the 
large-distance behavior 


^loss =r\A^\\ (36) 

where the recombination loss probability V equals (com¬ 
pare with Ref. [44) 1 


Fs{R) 


R^oo 


[A“(s) e-^^^ + A°'^\s)F^^] 



(29) 


V = {1- e-^’*) 


|1 - Sll^p ' 


(37) 


where a normalisation to unit flux is introduced for later 
convenience. For the real values of s, due to the repul¬ 
sive effective potential the wavefunction Fs{R) 

vanishes for i? —)• 0 (in the absence of three-body res¬ 
onance), and |A“(s)| = |A°'^*(s)|, i.e., the scattering is 
purely elastic. The losses thus come exclusively from the 
Efimovian sector. There, the strongly attractive effective 
potential —Sq/R^ gives rise to logarithmic waves at small 
hyperradius (i.e., in the limit where all three particles are 
close): 


The last ingredient is that the large-distance in¬ 
amplitude A™ is determined by the projection of 
the incoming plane-wave onto the Efimovian sector. 
More precisely, the hyperangular overlap (^isg 1')/'^°^) = 
V 2 f can be evaluated for i? —>■ 00 us¬ 

ing the stationary-phase method, with a result of the 
form 

j^in(O) ^ ^out(O) (38) 


(i?) ~ -k A°J^^ e-^^° 


m 

R — ^0 L V ^"^0 

(30) 

where a unit-flux normalisation is again introduced for 
later convenience. The three-body contact condition 
Eq. (18) becomes 




(31) 


then, Ag" —must vanish, because ((()s|'0~ has 
to behave like a purely outgoing wave at i? —)■ 00 , by 
definition of the scattering state ip [cf. Eq. (HU)]- This 
yields 

A“ = 2^/2 ^^/2g*5V4 ^-1/2 (3g) 


where 




( 32 ) 


^ The subscript 3 is used here since the subscript 2 is traditionally 
reserved for the atom weakly-bound-dimer channel |26ti28| . 
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after the hyperangular average: 
64 TT^ cos^(j) 


Inserting this into Eqs. (36) and (22), (j)s drops out 
■ the hyp 

K{k) = 


1-kiiP 


(40) 


fl — 6 - 

|1 — sii^p 

The final expression Eq. © follows immediately. 
B. Finite scattering lengths 


We turn to the general case where the scattering 
length(s) are not restricted to the unitary limit Eq. ([^. 
We will see that the final expression for the loss rate 
Eq.(§ remains valid, provided the definition of sn is ap¬ 
propriately generalized. Since the scale invariance of the 
two-body contact conditions is broken, the separability 
in hyperspherical coordinates no longer holds. Instead, 
we follow an S'-matrix approach. 

We consider the state 'I'l that physically corresponds 
to a stationary triatomic flow injected at the origin of the 
six-dimensional space (i.e., at R = 0) that gets partially 
reflected back and partially transmitted towards infin¬ 
ity. More precisely, is defined as the solution of the 
Schrodinger equation (16) with energy E = k'^/m sat¬ 
isfying the two-body contact condition and having 
the asymptotes 


4-1 (R) ~ 

fl-J-O 


{kRy^» +sii{kR)-^^° 
R 2 



’ 2 S 0 


m (41) 


and 


4'i(R) ~ S3ie 

R—^oo 


ikR 


m 1 
2kR W 


$3(0). (42) 


Together with the normalization ($ 3 [$ 3 ) = 1, this de¬ 
fines the reflection and transmission amplitudes sn and 
S31, as well as the function $ 3 (r 2 ) (up to multiplica¬ 
tion of S 31 and $3 by arbitrary phase factors e*''' and 
6 “*'-', respectively). Equation (41) is applicable in the 


scale-invariant region^ which we can define by i? <C 
mindoAel) 1/^) in case i and R ^ min(|aAB|) IobbI) 1/^) 
in case ii. Equation ( |42[ ) is applicable in the asymp¬ 
totic region: R 3> min(|aAB|) 1/^) in case i and R ^ 
min[max(|aAB |5 IobbI), 1 /fc] in case ii. 

Then, let {$„(r2)}^>^ be an arbitrary orthonormal set 
of functions such that {$n( 0)}^>3 forms an orthonormal 
basis [for the scalar product (.[.) defined in Eq. (27)]. A 
complete set of incoming and outgoing asymptotic states 
can be defined as 







2kR R 2 


(43) 


n > 3 , (44) 


and V-™* = (V^“)* for any n in the set C = {1} U {n; n > 
3}. The terms incoming and outgoing are meant with 


respect to the intermediate region contained in between 
the scale-invariant and asymptotic ones. 

For an arbitrary solution E of Eqs. (16) and dll 


in- and out-amplitudes can be defined by 


the 


fl-s-O 




A“ yy + A°“‘ 1 /:™* 

(45) 

n>3 

(46) 


The out- and in-amplitudes are linearly related: 

= E 


A° 


^nm-^rn 


(47) 


TTL^C 


where the matrix Snm is unitary and symmetric, as shown 
in Appendix [A| Furthermore, Snm is independent of the 
three-body and inelasticity parameters Rq and ij, and 
only depends on ka ab, fcuBB, and the mass ratio; indeed, 


we did not impose that E satisfies Eq. (18) as the three- 
body boundary condition. 

This problem of an a priori infinite number of coupled 
channels actually reduces to only two channels. Indeed, 
the channels 1 and 3 decouple from the other ones, i.e. 
Eq. (34) remains valid. To check this, first note that 
s„i = 0 for n > 4 by construction, cf. Eq. ( [4^ . Further¬ 
more, Sn 3 also vanishes for n > 4, because the state with 
a purely incoming wave in channel 3 (denoted by $3 in 
App .0 is a linear combination of d-i and $ 3 . 

The rest of the reasoning closely follows the a = 00 
case. Equations ( |31|32|33| hold, and hence also Eq. ( [^ . 
The s matrix being unitary, Eq. (36) follows. It re- 


3 to the projection of the incoming 


mains to relate A 
plane wave onto channel 3. For i? —>■ 00 , the overlap 
($ 5 I' 0 ^°^) can again be evaluated using the stationary- 
phase method, with a result of the form Eq. (38). On 


the other hand, the overlap (^sIV-) behaves at large R as 
{2kR) R~‘^ plus an outgoing wave. Since 
($ 5 I' 0 ~V'^°^) still has a purely outgoing behavior at large 


A, A™ — A™ = 0, which hnally gives Eq. (39) with 


replaced by $3. After hyperangular averaging, $3 drops 
out of the final expressions Eqs. (40) and (§. This hap¬ 
pens because we consider a non-degenerate gas at equi¬ 
librium, whose momentum distribution follows the Boltz- 
man law, so that the three-body momentum distribution 
depends only on the center-of-mass momentum and on k 
but not on k\ in general the functional form of $ 3 (^ 2 ) does 
play a role, see, for example, the study of non-equilibrium 
effects in Ref. HSl. 


C. Analogy with an interferometer 

It has been noted imiiH] that the loss peaks can be ex¬ 
plained by multiple reflections of the hyperradial wave off 
the intermediate region leading to a resonant denomina¬ 
tor under the integral in Eq. 0 - This behavior becomes 
transparent if we observe that the considered three-body 

























FIG. 3. A three-body wave arriving from large hyperradius 
R with amplitude ^ 3 " in the triatomic channel can follow 
various pathways before it either returns to large R, or gets 
lost at i? = 0 by turning into an atom and a deep dimer. 
There is a formal analogy with a Fabry-Perot interferometer 
with two partially reflecting mirrors. 


inelastic scattering problem is formally analogous to a 
simple interferometer with two partially reflecting mir¬ 
rors, see Fig.|^ The first mirror (located at intermediate 
R) has reflection and transmission amplitudes given by 
the 2 x 2 matrix 


f fill Sl3 
\S31 S 33 

The second mirror (located at i? = 0 within the zero- 
range model) has the reflection amplitude A which de¬ 
pends on the three-body parameter and the inelasticity 
parameter [cf. Eq. ( |32[ )]; transmission through this mir¬ 
ror corresponds to the three-body loss process, and hap¬ 
pens with probability 1 — e~^'^. 

It then becomes clear that the loss probability is mod¬ 
ulated by the interference between the different path¬ 
ways corresponding to multiple reflections by the two 
mirrors. More precisely, rewriting the term 1/(1 — snA) 
in Eq. (35) as X;„>o (sii^)", the n-th order term corre¬ 
sponds to the pathway with n reflections by each mirror. 


This is the origin of the term |l-|-(fci?o)~ | 


)-2iSog-277 gii|-2 _ 

|1/(1 — in the final expression Eq. ([^. This also 

clarifies why the dependence of the loss rate on Rq and 
r] is known analytically, Sn being independent of these 
parameters. 

The interferometer analogy also physically explains 
why the energy-dependent event rate constant K(k) has 
the upper bound 


K{k) < = 


64 TT^ cos^(f) 
mk'^ 


(48) 


This bound is a manifestation of the fact that the loss 
happens through a single Efimovian channel at short dis¬ 
tance (channel 1 ), which is coupled to a single large- 
distance channel (channel 3). The bound is reached when 
a perfect destructive interference leads to = 0. The 
lost flux (^loss is then equal to the incoming flux from in¬ 
finity |A™p (which is determined by a projection of the 
incoming plane wave, as we have seen). After thermal 
averaging, Eq. (48) implies 




32 TT^ cos^cj) 
m^{kBT)^ 


(49) 


where the inequality is now strict since K{k) = K-caax{k) 
cannot hold for all k (and we restored the h dependence). 

As a side remark we note that in the case 
where the B particles are fermionic, with me/mA S 
(13.60696 ...; 75.99449 ...) so that the Efimov effect oc¬ 
curs in the total angular momentum L = 1 subspace 
but not in higher L subspaces [Ml ESI 07], the bounds 
Eqs. (48) and (49) are multiplied by 3, due to the 3 
Efimovian sectors with angular momentum projections 
M = —1, 0 and -1-1. 

Let us now briefly address the low-energy limit. When 
fca —>■ 0 “, I Sill 1 , he., the first mirror becomes nearly 
perfectly reflecting. If 77 is small enough, the second mir¬ 
ror is also of good quality, and the finesse of the interfer¬ 
ometer becomes sufficiently high to observe resonances 
in the variation of the loss rate with scattering length. 
The resonances occur when a is such that there exists an 
Efimov trimer of vanishing energy. This happens when 
siiA « 1 , i.e. when the interfering pathways nearly lead 
to a divergence of (1 — SnA)”^ [IS]. This naturally leads 
to peaks in the loss rate constant K versus magnetic field, 
which remain visible and approximately unshifted after 
thermal averaging at low enough temperature. 

Finally, let us discuss the temperature-dependence of 
the loss rate at the unitary limit in case i (see also |44) j. 
From Eqs. ([^ and (10), one finds a qualitatively dif¬ 
ferent behavior depending on the order of magnitude of 
knl = When the mass ratio ms/mA is suffi¬ 

ciently large, as in the Li-Rb mixture discussed above, 
g-TTso jg gjgiaii compared to unity, so that can be ne¬ 
glected to a good approximation in Eq. ([^, leading to 
iF « (1 — e~'^^)K^ax oc 1/T^. In the language of the 
interferometer, this means that the first mirror is almost 
transparent, so that interference effects are negligible. 
The situation changes when ms/mA is not large, e. g. 
for the mixtures ®Li-^Li, ®®Rb-®^Rb, or '^°K-®^Rb (for 
which non-overlapping Feshbach resonances are indeed 
available [30])• Then, there is significant reflection from 
the first mirror, and interference effects lead to noticeable 
log-periodic oscillations of K as a function of T. To 
observe an entire period of these oscillations experimen¬ 
tally is challenging, since this would require changing T 
by the large factor However, some T-dependence 

of K may be detectable. 


D. Relation with previous work 


Let us comment on the similarities and differences be¬ 
tween the approach presented here (and introduced in 
Ref. [TU] in the BBB case) and previous work on three- 
boson scattering within the zero-range theory [261428) . 
Formally, it is possible to derive Eq. (1^ starting from 
the ^'-matrix formalism of Refs. j26H28| T However, the 
present approach is physically more transparent. A key 
point is that we directly constructed the relevant large- 
distance triatomic channel and its hyperangular wave- 
function $ 3 ( 17 ), cf. Eqs. ( 41|42 ). To this end, we consid- 
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ered the wavefunction d/i, corresponding to a triatomic 
flow injected at the origin. The idea of this wavefunc¬ 
tion d>i was already present in Ref. [22] where Efimov 
introduced the concept of s matrix connecting short- 
distance with large-distance hyperradial motion. How¬ 
ever Ref. |26| focused mainly on the case of negative 
total energy, where the only open channel at large dis¬ 
tances corresponds to motion of an atom relative to a 
weakly bound dimer, while triatomic motion is ener¬ 
getically forbidden. Braaten and Hammer [27] added a 
triatomic large-distance channel charaterized by a wave 
function independent of the hyperangles. The resulting 
s matrix provides a suitable framework for studying the 
three-body recombination in a Bose gas with finite a in 
the zero-temperature limit. Then, in Ref. |28j this for¬ 
malism was generalized to finite temperatures by using 
a complete set of long-distance channels with hyperan- 
gular wavefunctions defined by hyperspherical harmon¬ 
ics. A conceptual difficulty with this construction is that 
these large-distance channels only decouple for |a|, 
so that one has to formally work with a finite a. In 
contrast, one would expect physically a smooth depen¬ 
dence on 1/a in the interval (—oo;0]. This expectation 
is confirmed by the present construction. The shape of 
the function $ 3 (r 2 ) depends only on ka and interpolates 
between two limits: for small k\a\ it is a constant inde¬ 
pendent of Cl and for large k\a\ it tends to (j)isg{^), the 
asymptotic region in this case being defined by i? 3> 1 /fc. 
Incidentally, this illustrates the breakdown of the adia¬ 
batic hyperspherical approximation for k\a\ ~ 1 . 


IV. CALCULATION OF THE FUNCTION sn 


/bb defined by 


/ab(x) = 47 r lim 

y-fO 

/bb(X) = 47 r lim y vki, (50) 

y —>0 

where (x, y) are the Jacobi coordinates given in Eqs. ( [ll] ), 
while (X, Y) is another set of Jacobi coordinates defined 
by 


cos 6 » X = r 2 - (ri -f r 3 )/ 2 , 


(51) 


2sin0 Y = r 3 — ri. 

The STM equations conserve the angular momentum 
and, in the considered regime where the mass ratio is not 
very large, the Efimovian solution appears only in the 
channel with zero angular momentum, so that /ij(x) = 
fij{x). For this reason the contribution of higher angular 
momentum channels to the loss vanishes in the zero-range 
limit, in contrast to the three-body recombination into a 
weakly-bound state for positive a when ka is not small 
pH] . For more details on the derivation of the STM equa¬ 
tions in the general case of different masses and scatter¬ 
ing lengths see, for example. Ref. [ID]. The asymptotic 
behavior of 'I'l at small hyperradius given in Eq. (41) 
translates into 


fAB{x) 


f Cab^ ikxy‘'° + sii{kx) 

>0 \Cbb J X 


(52) 


where Cab and Cbb are numerical coefficients given in 
App. I At large x both functions fAB{x) and /BB(a;) 
should represent outgoing waves: 


/AB(a:)l ^ expjikx) 

fBB{x) x—>-oo 


(53) 


Our computational method for the universal function 
sii(fcaAB, fcoBB) consists of finding the three-body wave- 
function 'hi defined by Eqs. (41) and (42) numerically. 
We consider case u, while case i is treated similarly mod¬ 
ulo the simplifications listed at the end of this section. 
Analytic results for sn near the unitary limit are pre¬ 
sented in Appendix [P] 


The Schrodinger equation together with the two-body 
contact condition for the wave function 'I'i(R) can be 
reduced to a set of coupled integral Skorniakov and Ter- 
Martirosian (STM) equations for the functions /ab and 


We mention that at the unitary limit, the solution is 
known analytically and fij are expressed in terms of the 
outgoing Hankel function 


fAB{x) 

fBB{x) 


/Cab j 


3 / 


gSQTl 


(54) 

which, matched with Eq. (521, gives the limiting expres¬ 
sion Eq. (10). 

In order to calculate sn for arbitrary {kaAB,kaBB) 
we switch to momentum representation where the STM 
equations write 


{\/P^ — P — f0+ — l/aAB)/AB — -bfc2,0/AB — ifc2,e/BB = 0 

,_ .. (55) 

{\^P^ — — f0+ — 2 sin0/aBB)/BB — 21/^2 ^/ab = 0, 

where fijip) = J ftj{x) ex.p{-ipx)d^x, 

1 f°° ')^' \ f P'^ “^PP'P^ ~ ^ ~ 

TT sin 2a Jq p \p''^ — 2pp' sin a + p'^ — k'^ cos^ a — i0+ J ’ 


Lk^afip) = 


(56) 
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(j) and 6 are the mass angles defined in Eqs. ([^ and (|^, and the inclusion of a small positive 0+ indicates that k is 
slightly shifted into the upper complex halfplane thus fixing the branch cuts of the logarithm and square root and 
ensuring the presence of only the outgoing wave (ISSl) in the solution. The boundary condition (|52|) now reads 


/ab(p) 
/bb (p) 


f Cab\ 27r^so 


p-)- 


OO yCBB/ sinh(7rso/2) 


1 


{p/k)- 


r(i-zso) 


+ Sll 


1 


{p/ky 


r(i + *so) 


(57) 


We then perform a complex scaling transformation. We use the fact that sn depends only on the products oab^ 
and obb^ and replace in Eq. (55) fc —>■ i, Oy —>■ —iaijk. Equivalently, this can be done by rotating the integration 


contour over p' in Eq. (56) to the negative imaginary axis. The obtained STM equations, 


{'Jp^ + 1 — i/ aAB^)/AB — ^-i,<^/ab — i-ye/BB = 0, 
(^^2 2ism0/aBBk)fBB - 2T-i,e/AB = 0, 


(58) 


do not contain singularities in the kernels and can be efficiently solved numerically. We deduce sn from matching the 
solution with the asymptotic form 


7ab(p)' 


fCAB\ 

_/bb(p)_ 

(X 

p—>-oo 

^^Cbb j 




r(l + iso) 


(59) 


for any given set oab^ and obb^- More precisely, we 
use the following procedure. We introduce new un¬ 
known functions 5ab(p) and 5bb(p) such that fijip) oc 
[p“2-*so _|_ gj;(p)] [the proportionality coefficients follow 
from Eq. (57) but are irrelevant for the determination of 
sii] and require that for large momenta gij{p) oc 
Equations (58) become a set of linear inhomogeneous 
equations for qab and gsB, which we solve by discretiz¬ 
ing p and inverting the corresponding discrete analog of 
the operator in the left hand side of Eqs. (58). The 
final result is given by sn = e“’^®‘’r(l -|- i^/r(l — 
iso) linOp-j-oo (p)p^~*®“ independent of the choice gxB 
or pbb- 


Case i is obtained by setting /bb = Cbb = obb = 0 
in the above analysis and omitting the second lines in 
the STM equations (551 and (58). Various properties of 


the ABB system in case i have been studied by using 
the STM equation derived in this case from the effective 
field theory . The STM equations in both cases i and 
ii have been employed in a recent study of the single¬ 
particle momentum distribution of heteronuclear Efimov 
trimers m- 

Another point which is useful to mention is that 
Braaten et al. [28] have developed a slightly different 
method for calculating sn in the case of three identi¬ 
cal bosons. They also use the STM equation but deduce 
sii from solving the atom-dimer scattering problem for 
positive a above the dimer break-up threshold. We notice 
that their approach works well for small ka > 0 whereas 
ours is optimal for large ka, either positive or negative. 
In fact, sii(-aABfc,-OBBfc) = /sl^iaABk, aBBk). 


This relation follows from Eq. (59) and the observation 
that if fij is the solution of Eqs. (58) for a^, then is 


the solution for ay = —ay. 


V. SUMMARY AND OUTLOOK 


In this paper we have obtained the finite-temperature 
three-body loss rate constant within the zero-range the¬ 
ory for the ABB system, with uab < 0, and no BB in¬ 
teraction (case i) or obb < 0 (case ii). For a given mass 
ratio, we expressed the rate constant in terms of tempera¬ 
ture T, three-body parameter i?o, inelasticity parameter 
rj, and a universal function sn that depends on ka^Bi 
as well as on fcoBB in case ii, where k is the three-body 
collision momentum. We developed a numerical method 
based on complex scaling for computing the function sn 
and perform this calculation for two experimentally rel¬ 
evant cases: ®Li-^^^Cs-^^^Cs and ®Li-®'^Rb-®^Rb. The 
knowledge of sn reduces the problem of computing the 
loss rate to a simple thermal averaging integral over k for 
any desired T, Rq and g. We expect that these results, 
combined with experimental data, will be useful for pre¬ 
cise tests of universality and determinations of Rq and 
g. For ®Li-^^^Cs-^^^Cs we find that inclusion of the CsCs 
interaction leads to a significant enhancement of the loss 
rate constant. This is in spite of the fact that the scaling 
factors in cases i and ii are very close for this large mass 
ratio. The enhanced loss rate is likely to be explained 
by an enhanced probability to find two Cs atoms close to 
each other in the CsCs-interacting case. 

In deriving the loss rate coefficient we explicitly con¬ 
structed the hyperangular wave function $ 3 (r 2 ) corre¬ 
sponding to the long-distance three-atom channel which 
is connected by a unitary two-by-two matrix with the 
Efimovian wave at small hyperradii. The corresponding 
S'-matrix formalism smoothly connects with the exactly 
solvable unitarity limit. The three-body loss problem re¬ 
duces to a simple Fabry-Perot interferometer with two 
partially reflecting mirrors. 
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Three-body systems at zero total energy are known 
to be analytically solvable in the zero-range approxima¬ 
tion if there is only one relevant scattering length (three 
identical bosons gOlllHHSI], ABB system in case i [33], 
ABB system with fermionic B particles jH]). A theoret¬ 
ical challenge that we can formulate now is to generalize 
these approaches to the ABB system in case ii at vanish¬ 
ing total energy. In particular, it would be interesting to 
obtain an analytic expression for the loss rate constant 
at zero temperature (the dash-dotted line in Fig. [^. 
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Appendix A: Basic properties of the s matrix 


Here we derive some basic properties of the s matrix 
that were stated and used in the main text. 

Let us first give a proper definition of the s matrix. 
In addition to the state d*! defined above, let us de¬ 
fine the state for m > 3 as the solution of the 
Schrodinger equation with the two-body contact condi¬ 
tion Eqs. whose asymptotic behavior contains an 

incoming wave of unit amplitude in channel m and a 
purely outgoing wave in the other channels; the coef¬ 
ficients of the outgoing waves then define the column 
(snm)nec of the s matrix: 


R^oo ^^ 

n>3 


(Al) 


A useful lemma is that for an arbitrary state F solving 
Schrodinger’s equation with the two-body contact condi¬ 
tion Eqs. ([I6|T7| ), the in- and out-amplitudes [defined in 
Eqs. (45|46)j are constrained by 


= (A2) 

n^C nGC 


This follows from the conservation of probability. More 
precisely, the flux through a hypersurface S [defined as 


in Eqs. ( 19|20 |] is independent of S. Taking for S a 
hypersphere of very small or very large radius, the flux 
respectively equals or — 

\K?)- 

A useful co nseque nce is that if and F' are two solu¬ 
tions of Eqs. ( 16|17 1 with the same in-amplitudes, then 
their out-amplitudes are also equal. This follows from 
applying the above lemma toF — F'. 

T he l inear relation between out- and in-amplitudes 
Eq. (471 then follows by noting that F and 
have equal in-amplitudes and hence equal out- 
amplitudes. 

Let us now check that the matrix Snm is unitary. Phys¬ 
ically this comes again from the conservation of probabil¬ 
ity. More precisely, applying the lemma Eq. (A2) to the 


state F — -|- a where n ^ n' and a is an arbitrary 

complex number, we get 


1 + lap = P" 

mGC 


■ OL Sri 


(A3) 


Taking a = 0 yields 


= (A4) 


Thus (A3) simplifies to B,e{aJ2m = 0 for any 

a, which implies 


'y ^ ^mn^rnn' 0- (AS) 

m 

Let us now check that Snm is symmetric. Physically, 
this follows from time-reversal invariance. More pre¬ 
cisely, we consider the state 4'*, and notice that it has 
the same in-amplitudes than X(mGC ^mn hence also 
the same out-ampltiudes. This yields s • s* = 1, and 
hence, since s is unitary, s = s"^. 


Appendix B: Determination of so, Gab and Gbb 


The STM equations can be used to calculate sq and 
the ratio of Cab and Cbb- To this end, it suffices to 
consider the scale-invariant case where the interaction is 
at the unitary limit and the energy is zero. 

Let us first consider case ii. A solution of the STM 
equations (551, for fc = 0, is given by the ansatz m 


■/ab(x)1 ^ /Cab\ ^ 
/bB^x) \Cbb/ 


Indeed, we have /p^ = A(so, where the 

function A is defined by Eq. (§. The STM equations then 
become a 2x2 homogeneous system of linear algebraic 
equations for Cab and Cbb: 


(l-\{sQ,(j)) -A(so,^')\ (Cxb\ 
\-2X{so,e) 1 Jl^CBByl" 


(Bl) 
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The requirement that the determinant of the 2x2 ma¬ 
trix in Eq. ( |B1[ ) vanish gives the implicit equation for sq 
E q. 0. The ratio between Cab and Cbb is then fixed by 
Cbb = 2A(so, 0)Cab and we can choose these coefficients 
to be real. 

In case i, one should formally set Cbb = 0 in the above 
analysis. 

Finally we note that, even though all that is used for 
the computation of sn in Sec. |TV] is the ratio of Cab 
and Cbb, their absolute values can also be determined, 
thanks to the expression of the normalized wavefunction 
given in Appendix [C| Indeed, Eqs. ( 41|50|C3|C4 ) 
yield straightforwardly 

= 47rsinh(^so|) > (B2) 


\Cbb J 


2so V^bb/ ’ 


and the hyperangular part writes 
</.,«„(J^) = (1 + P)Cab—^^^+Cbb—(C3) 


sin a cos a 


sin /? cos /? 


where 


ip{a) = sinh Sq 


a and j3 being hyperangles defined by 


X = R cos a 
y = R sin a 


and 


(C4) 


the absolute values of Cab and Cbb being determined by 
Eq. ((C7|. 


X = R cos /3 
Y = i? sin P, 


Appendix C: Normalised wavefunction and contact 
parameters of a heteronuclear Efimov trimer 

In this Appendix, we consider the three-body bound 
states, i.e. the negative-energy solutions of the zero- 
range model [defined by the three-body Schrddinger 
equation Eq. (16), the two-body contact condition 
Eq. ([^, and the three-body contact condition Eq. 
at the unitary limit Eq. (|^ and without losses {rj = 0). 
We provide the analytical expressions of their normalised 
wavefunction and their contact parameters. We consider 
case ii throughout this Appendix. The content of this 
Appendix also applies to case i provided Cbb is formally 
set to zero. 


the Jacobi coordinates {x,y) and {X,Y) being defined in 
Eqs. 0 and ( [^ . 

Imposing the two-body contact conditions Eqs. 0 
yield s th at Cab and Cbb satisfy the same matrix equation 
Eq. (B1) than Cab and Cbb, which again yields Eq. 0 
for So, as well as 


Cbb = 2A(so, 0)Cab- 


(C5) 


The spectrum is 


En = -?_4e2"gr(i+iso)Aog-n2V-o^ ^ (C6) 

tuRq 


1. Wavefunction 


The trimer’s wavefunction writes 1241 




(Cl) 


where the hyperradial part is proportionnal to a Bessel 
function 


f{r)=j^k,so{V^\R) 


as obtained by imposing the three-body contact condition 
Eq. (pi). 


The normalization of the wavefunction can be done 
analytically, generalizing |53j to the heteronuclear case. 
For the hyperangular part, we take {4)iso\4'iso) = 1 for 
the hyperangular scalar product introduced in Eq. (27). 
This leads to 

(2Cab Y C^'q)Q{'x/2) + 2CabQ{^) + 4CabCbbQ(^*) = 1, 


(C2) where 


(C7) 


Q[<^) = - ^cosh fsoTr) sinh(so a) — 2a sinh fso —) cosh(so a) 

sn sm 2a L \ 2) \ 2 J 


So sin 2a . 


(C8) 


r 


and 


Q{tt/2) = lim Q{a) = 47r^[sinh(so7r)/so — tt] 

ol ^ t : 12 


In order to normalize the hyperradial wavefunction, 
it is convenient to consider that the three particles 
are subject to an external harmonic trapping poten- 
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tial of vanishing frequency |53j . so that one can impose 
/ <Pri(Pr 2 (Pr 3 \^{Yi,V 2 ,vz)\^ = 1, where 'I' equals ^/’(R) 
times a center-of-mass wavefunction 'i/'Cm(C) norma lized 
to f d^C'|V'CM(C)p = 1. Due to the Jacobian Eq. (23), 
this gives dRR |F(i?)p = cos“^ (j). The integral over 
R has a known expression [54) . which finally yields 


N = 


12m\E\ sinh(so7J‘) 

SqT^ COS^(j) 


(C9) 


2. Contact parameters 


From this normalized wavefunction, it is straightfor¬ 
ward to deduce the contact parameters of an Efimov 
trimer, i. e. the partial derivatives of its energy with re¬ 
spect to the scattering lengths taken at fixed three-body 
parameter. Indeed, as shown in Ref. [^, (dE/daij)^^ is 
proportional to the norm of the function Aij appearing 


in the 2-body contact-condition Eq. (17). This yields 


dE 


9(-l/aAB) 

dE 


i9 (—I/obb) 


2 ^ a 3 ta'iih(so7r) sinh^ (so7 ’‘/2) 


— ^AB 1/ — l^TT' 


So 


-n 2 \E\ 3 tanh(so7r) sinh (so7r/2) . 

— Cbb \/ - iOTT - Slllc/. 

\ m So 


(CIO) 

(Cll) 


This generalizes to the heteronuclear case the result obtained in [SSI for three identical bosons. 


Appendix D: Expansion of sn around the unitary limit 


This Appendix concerns leading-order corrections to sn near the unitary limit. We start by stating the results. In 
case ii, when /cIoabI and fc|aBB| tend to oo, we have the expansion 


Sll 


1 — IGtt^ tanh(so7r) sinh^^SQ—^ 




AB 


Cgg sin0 


yOAB^ Obb^ 


(Dl) 


r 


where is given by Eq. (10) a nd the coefficients Cab 
and Cbb are derived from Eqs. ( C5|C7 ). In case *, the 
expansion for fc|oABl oo is given by formally setting 


Cbb = 0 in Eqs. (D1|C7). 


We present the derivation in case ii (case i is treated 
similarly with obvious simplifications). We employ the 
same complex scaling procedure than in Sec. |IV[ now at 
the level of the Schrodinger equation. Namely, we assume 
that the function Z i—4'j(^R,)^ defined a priori for real 
positive Z, can be analytically continued to the quadrant 
0 < ArgZ < 7r/2, and we consider the scaled wavefunc¬ 
tion ^i(R) = 4'i(iR/A:) where R has real coordinates. 
The Schrodinger equation —ArJ^i = k^'^i then be¬ 
comes —= —^ 1 . Furthermore, the large-distance 
behavor oc gives Ti oc e~^ after analytic contin¬ 
uation. Finally, because 'I'j^ satisfies the two-body con¬ 
tact conditions with scattering lengths (oabjObb), J'l 
satisfies the two-body contact conditions with the scaled 
scattering lengths 


(oab, Qbb) = (qab, obb) k/i. 


(D2) 


The short-distance behavor dti oc {kRy^° -|- sn (fci?)“®®° 
turns into Ji*! oc which we rewrite as 

oc {R/Ro)~'^^° — (i?/i?o)*®“ where Rq can be viewed as 
an effective (in general complex) three-body parameter. 


scaled wavefunction corresponds to a bound trimer state, 
of fixed energy — 1 /m, with imaginary scattering lengths 
and complex three-body parameter. 

At the un itar y limit, the rescaling of the scattering 
lengths Eq. (D2) has no effect, and J'l is the wavefunc¬ 
tion of a standard Efimov trimer. Hence its effective 
three-body parameter Rp is real, and is related to its en¬ 
ergy, — Ijm, as in Eq. (C6). This allows to retrieve the 
expression Eq. (10) of s“. 

If we move slightly away from the unitary limit by turn¬ 
ing on small finite inverse scattering lengths, Rp has to 
shift in such a way that the scaled energy remains fixed. 
But the partial derivatives of the energy with respect to 
the inverse scattering lengths are given by the contact 
parameters, computed in App. 0 The partial derivative 
of the energy with respect to the three-body parameter, 
on the other hand, is obtained easily from the simple re¬ 
lation between energy and three-body parameter valid at 
the scale-invariant point. This yields the desired result 
Eq. dml). 


related to sn by sn = —Rp 


Summarizing, the 


Appendix E: Three identical bosons 

The case of recombination between three bosons in the 
same internal state, which was the subject of Ref. m, 
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can be recovered from the present article, modulo mi¬ 
nor modifications given in this Appendix. One sets 
iriA = m-B = m so that ij) = 9 = tt/Q, and oab = obb =: ct 
(case i does not exist any more). In Eq. §) the right 
hand side is multiplied by 2. For the symmetrization 
and normalisation of the plane wave, in Eq. ( |14[ ), the 
term {1 + P) 1^/2 is replaced by Pa) /'s/s!, where the 
sum over a now runs over the 3! permutations of the 
three particles. As a consequence, the right hand side 
of Eq. (391 (and of its finite-a generalization) gets mul¬ 
tiplied by v^S. Furthermore, in Eq. (22) the right hand 


side contains an additional factor 1/3. The additional 
factors cancel out in the final result for K and Eq. © is 
unchanged. 

As far as the calculation of the function sn is con¬ 
cerned, Sec. |IV| remains entirely valid in the case BBB, 
but gets simplified by observing that /ab = /bb and 
Cab = Cbb which leads to the usual single STM equa¬ 
tion. Similarly, all analytic results in Apps. [gand[g 
reproduce the known BBB ones (in this case Cab = Cbb , 
and di/aAB + '^^i/obb should be replaced by di/a)- For 
completeness we also mention the relation between a_ 
and i?o in this case, a_ « — 1.017i?o. 
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